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Abstract 

We employ notions familiar from supersymmetry for constructing the one-loop func- 
tional of general quantum field theories with bosons and fermions (spin < 1/2). To 
demonstrate the advantages of such an approach for calculating one-loop divergences, 
we analyse a simple Yukawa theory with two different versions of a super-heat-kernel 
expansion. These methods also simplify the calculation of one-loop divergences for non- 
renormalizable meson-baryon Lagrangians occurring in chiral perturbation theory. 
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1. The one-loop functional of a general bosonic or fermionic theory can be expressed in 
terms of the determinant of a second-order differential operator. A common procedure in 
theories with both bosonic and fermionic degrees of freedom is to first integrate out the 
fermions and then treat the resulting bosonic theory. 

The purpose of this note is to demonstrate the usefulness of treating bosons and 
fermions on the same footing. Although this has been common practice in supersymmetric 
theories from the early days of SUSY, we emphasize the advantages of such an approach 
also for non-supersymmetric quantum field theories with arbitrary numbers of bosonic 
and fermionic fields. Here, we consider only fields with spin < 1/2. 

The method turns out to be especially useful for calculating the divergent parts of one- 
loop functionals occurring in the renormalization program of chiral perturbation theory 
with mesons and baryons |l|, @, |3|, ||. In the standard approach, the bosonic loop, the 
fermionic loop and the mixed loop (boson and fermion lines in the loop) are treated 
separately |5], |(J. This requires a cumbersome investigation of the singular behaviour of 
products of propagators, because the mixed loop does not have the form of a determinant, 
as the purely bosonic or fermionic loops. The SUSY-inspired treatment, on the other 
hand, reduces the problem to simple matrix manipulations, in complete analogy to the 
familiar heat-kernel expansion technique for bosonic or fermionic loops. 

We demonstrate the simplicity of this approach in the case of a Yukawa theory, where 
the divergences can relatively easily be extracted by standard Feynman diagram tech- 
niques. With applications to baryon chiral perturbation theory in mind 0], we will actu- 
ally analyse the Yukawa theory in two different ways. In the first method, we square the 
fermionic differential operator in the usual way to arrive at the superspace version of an 
elliptic second-order differential operator (in Euclidean space). In the second method, we 
set up the heat-kernel expansion directly for the original supermatrix. Since we are not 
aware of a discussion of this technique in the literature, we will be more explicit for the 
second method. 

2. We start from a general Euclidean action 



for n B real scalar fields <pi and rip spin 1/2 fields ip a . To construct the generating func- 
tional Z of connected Green functions, we couple these fields to external sources ji 
(i = 1, . . . ,n B ), p a , p a (a = 1, . . . ,n F ), 




(1) 




(2) 



where we have used the notation 



fif + ipp + pip : 



/ 



d d x (ji<fi + ip a p a + p a ip a ) . 



(3) 
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The normalization of the functional integral is determined by the condition Z[0, 0, 0] = 0. 
We denote by <p c \, ip c \ the solutions of the classical equations of motion 

5S 5S 5S 

difi bip a dip a 

With fluctuation fields £, r\ defined by 

<fi = <£cl,i + & 

Ipa = 1pcl,a + Va , (5) 

the integrand in (|2|) is expanded in terms of £, rj, rj. The resulting loop expansion of the 
generating functional 

z = z L=0 + Z L= i + . . . 

starts with the classical action in the presence of external sources: 

Z L=0 = S[(p ch ip ch ip cl ] - j T (p d - 4> cl p - pip c i , (6) 

where the classical fields are fixed by the external sources through (||) . The one- loop term 
Zl=i is given by a Gaussian functional integral 

e~ z ^ = J [d^dr]dr]} e -s w lv*,1>*$*£,vif>\ ; (7) 

where 

S {2) [(fi ch ip ch ip d ; 6, V,v} = / d d x £ (2) (ip ch ip ch ^ cl ; £, 77, 77) (8) 



is quadratic in the fluctuation variables. Employing the notation introduced in (Q), 
takes the general form 

= ~(fAZ + rjB V -r) T B T ri T + fTr] - V T T T ^ + m ~ f^f) , (9) 

where A, B, V are operators in the respective spaces; A = A T and B are bosonic differen- 
tial operators, whereas V and V are fermionic (Grassmann) operators. They all depend 
on the classical solutions (p c \, ip c \. The second explicitly symmetric form of in (|9|) will 
be used for the definition of the supermatrix operator K in (|TT1). 



The standard procedure for the evaluation of (|7|) is to integrate first over the fermion 
fields 77, rj to yield the bosonic functional integral 

Z L=1 — r|p+ R / —h(; T (A—TB~ 1 T+r T B~ 1T T T )!l 



e -z L=1 =detB J ^] e - 
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In this way, we obtain the familiar result 
1 



J L=\ 



lndet(A-r£- 1 r + r T B- 1T r T ) -In det A) - (In det B - In det B 



- Tr In 4- ~ Tr In — + - Tr ln(l - A^TB^T 
2 A S 2 1 



- Tr In — Tr In — Tr 

2 A B ^2n 



A 



-i 



A-^B-^V 7 ) 



YB^Y - A^Y 7 B- lT Y T ) n , (10) 



where 



An 



j=p=p=0, 



B :=B 



j=p=p=o 



denote the free-field limit of A and B, respectively. Recalling that are the 

scalar and fermion matrix propagators in the presence of external sources, the one-loop 
functional Zl = i is seen to be a sum of the bosonic one-loop functional | Tr ln(A/A ), 
the fermion-loop functional — Tr ln(B/B ) and a mixed one-loop functional where scalar 
and fermion propagators alternate. In order to determine the ultraviolet divergences that 
occur in the last term in (|T0|), the calculational inconveniences mentioned in paragraph 
1 are encountered. 

In the remainder of this letter, we discuss a procedure that allows us to identify the 
singular pieces in Z^ = i more directly. Inspired by supersymmetric quantum field theories, 
we reorganize the three parts of Zi, = \ into a more compact form, using the notion of 
supermatrices, supertraces, etc. (cf., e.g., Refs. f8|, |9j, |T0|1 ) . 

Assembling the bosonic and fermionic fluctuation fields into a multicomponent field 



in (H) can be written as 



SCO 



-\ T K A 



/ 


A 


Y 


-Y 




— T 




-B 




-Y 





V 


Y 


B 






where K is the supermatrix operator 

K 

For a general supermatrix of the form 

M = 



where a,b (a,/3) are bosonic (fermionic) variables, one defines the supertrace str M and 
the Berezinian sdet M as @, 




str M 
sdet M 



tr a — tr b , 

det (a — ab~ l (3)/ det b . 
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These definitions give rise to the relation 

sdet M = exp (str InM) , (12) 

in analogy to the one for ordinary matrices. 
Comparing this with Eq. QIDD, we obtain 

1 K 1 K' 

Z L=1 = -Str In — = -Str In— (13) 
Z K Z K Q 

with _ 

K' = \ y/Jl T fiB . (14) 

\^r T v bT / 

With our notation 

Str O = J d d x str(x\0\x) 

we are distinguishing supertraces with and without integration over Euclidean space. 
In ( |T^D we have introduced a mass parameter /i that guarantees equal dimensions for 
all entries in K' ([K f ] = [A] = 2). Although this quantity does, of course, not appear 
in any final result, it turns out to be quite helpful for the inspection of expressions at 
intermediate stages of calculations. 

It is seen that the three types of one-loop functionals (bosonic, fermionic and mixed) 
in Eq. (|10D are not distinguished any more in the representation (|13"D, where Zl=i is 
expressed in terms of a Berezinian. Generalizing the heat-kernel expansion to supermatrix 
operators, the determination of the divergent part of Z^ = i reduces to matrix algebra, in 
complete analogy to the purely bosonic or fermionic loop functionals. 



3. As long as we are only interested in those parts of the one- loop functional that are 
at most bilinear in fermion fields, we can reduce the supermatrix K' to the simpler form 



K"=[ /rT~~ -p ) , (15) 



izp r fiB 

such that the one-loop functional can be written as 

^ 1 = i strln |_i^ ln | + ... (16) 

The terms omitted are at least quartic in the fermion fields. 

The representation ( |i~6|) can be written in a more compact form by "squaring" the 
fermionic differential operator B. We first note that 

B B B B 

Tr In — = lndet — = lndet 75 — 75 = Tr In 75 — 75 = —Str In 
fi fi fi 
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Suppressing terms of higher than second order in the fermion fields in the rest of this 
paper, we get 

Z L -i = - Str In — - , (17a) 

with 

a=( A v^r 75 5 T5 \ 



In many cases (as the Yukawa theory to be considered below), this procedure brings the 
supermatrix differential operator A into the form 

A = —D^Dp + Y , (18) 

where = + and X^, Y are matrices in superspace. Generalizing the ordinary 
heat-kernel expansion to superspace, we will extract the second Seeley-DeWitt coefficient 
of the generic operator ( |I8[) to obtain the one-loop divergences including the terms bilinear 
in fermion fields. 

In order to avoid the pitfalls with 75 in d dimensions, the steps leading from Eq. (|T6|) 
to (|18D may be carried through by using an intermediate regularization that allows the 
evaluation of determinants in four dimensions. An example of such a regularization is 



provided by Eq. (IS), where the integration over the parameter r may be cut off at the 
lower end. In the following, we assume that X^, Y are independent of 75, as is the case in 
the Yukawa theory, and we return to dimensional regularization for ease of comparison 
with the standard methods to evaluate the loop integrals. In cases where 75 pertains (as 
in pion-nucleon effective theories), one may stick to the regularization just mentioned. 
In the proper-time formulation, the one-loop functional assumes the form 

= -if -Str (e---e-^ 



2 Jo T 

1 POO P 

= —— / — / d d x str (x|e" TA -e- TAo |x) . (19) 

2 JO T J 

To extract the divergences, we need the coefficient a,2(x,x) in the heat-kernel expansion 

00 

( x \ e - rA \x) = (Anr)~ d / 2 r n a n (x, x) . (20) 

n=0 

The divergent part Z^l x of the one-loop functional is then given by 



Z * Zl ~ (47r) 2 (d-4) ld * XStT 



a,2(x, x) — a?,(x, x) . (21) 



The derivation of the diagonal super-Seeley-DeWitt coefficient a 2 (x,x) is completely 
analogous to the ordinary case. Employing the method of Ball ||11||, we write 



, d d h 

i!(J„L-tA|il\/)J™\ / -ikx-rAikx 



(x\e- T *\x) = I d d k(x\e- T *\k)(k\x) = J J^y e e e 

d<lk e -r(fc 2 -2ifc M D M +A) 1 _ ^2) 



(27T) 
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With the dimensionless integration variable I = y^rk, we obtain from (p0|) 

d d l , 



a 2 (x,x) 



7T" 



^A 2 + ^(A/^A, + D tl AD v + D„A,A) 



and therefore, after integration over /, 

I i 

a 2 (x,x) = —X^X^ + -Y 2 



--DY - -YD + -D U YD U 1 , 
6 6 3 



(23) 



(24) 



with 



Xftv — [Dp, D v \ — dpX v — d u Xp + [Xp, X v 



By partial integration, the last term does not contribute to the divergence functional (|2T| ) 
and we arrive at the final result 



J d 4 x str a 2 (x,x) = J d 4 x str ^— X^X^ + ^Y' 



(25) 



4. We illustrate the formalism with a specific case, the Yukawa theory of one scalar and 
one fermion field. We consider the following Lagrangian in Euclidean space: 

£ = \v(-d 2 + M 2 )ip + + V(7A +m- gif)ip - j<p -ipp-p^p. (26) 

The one-loop divergences of this theory are due to the Feynman diagrams in Fig. |l|: 
the first two diagrams are contained in the bosonic one-loop functional, the next four 
are included in the fermionic counterpart, whereas the last two belong to the mixed 
functional, with both bosons and fermions running in the loop. 
The equations of motion are 



:-d 2 + M 2 )p cl 



■ A 3 T , 

J - gj Pel + 9W C \W C \ 



( $ + m)V>ci = p + gipdipd . 
In terms of the fluctuation fields defined in Eq. (|^), the fluctuation Lagrangian is 

£ (2) = + M 2 )i + r](0 + m)r l -g ^ clV + # cl ) , 



(27) 



M* = M 2 + K 2 cl , 



m = m- g(p c i 



Thus, the three entries of the supermatrix K" in Eq. (|15|) are 



A 



-d 2 + M 2 , B 



m 



-g V>ci • 
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Figure 1: One- loop Feynman diagrams with divergent parts for the Yukawa Lagrangian 
(|2"ED. Bosons (fermions) are denoted by dashed (full) lines. 



The supermatrix A in Eq. ( [[ 7b|) has the structure ([18]) with 



X, = I ° ^ ) , Y = f ^ ^ ^ ! ■ (28) 



In our simple example, without vector fields and without derivative interactions, the 
curvature X^ v does not contribute because X^X U = 0. The divergence Lagrangian takes 
the form 

2(47r) 2 (d-4) 

reducing the whole problem to a simple matrix multiplication. The final result in Eu- 
clidean space is 

£ dlv = {Aii)2 \ d _ 4) [\{M* - M 4 ) - 2(m 4 - m 4 ) - 2g%<p^<p d 

+ 9 2 M- ^ + 2m)^ cl } . (29) 

An explicit calculation of the diagrams in Fig. [1] reproduces exactly these divergences. 
The same result can be obtained with the background field method |12| . 



5. We now present an alternative method for calculating the one-loop functional by 
applying the super-heat-kernel expansion directly to the operator K' in (|I~4"D or to K" 



8 



in dig). For the Yukawa theory we are considering here, the previous method is simple 
enough. However, for non-renormalizable Lagrangians with derivative couplings occurring 
in chiral perturbation theory, the expressions can become more compact in the following 
approach. In addition, we wish to demonstrate with an explicit example that the super- 
heat-kernel expansion is perfectly well defined also in the linear version. 

We confine ourselves to the mixed functional in this case, i.e. to the part bilinear in 
the external fermion fields. Going through essentially the same steps as in paragraph 3 
and specializing immediately to the Yukawa Lagrangian (p6|), we obtain 



with 



V = -diag ((/ - *v/7«9) 2 + rM 2 , ii[i^/ + r(^ + m)}) , 

w = ^ TS ( I ^ l ) ■ (31) 



The appropriate decomposition of the exponential in (|30| ) can be performed by using 
Feynman's "disentangling" theorem ||13||: 



exp(V + W) = exp V P s exp ds W{s) (32) 

Jo 



with 



Vrxr„sV 



W(s) := e~ sV We 



and 

"1 00 rl fsi 



ri — "x, r 1 r s i rs n -i 

P s exp / dsW(s) := V / ds x / ds 2 . . . ds n W( Sl )W(s 2 ) . . .W(s n ) . (33) 
Jo n^cr J° J° 

Since we are only interested here in terms with the structure ip cl . . . ip c \, we pick out the 
part bilinear in W: 

str (e v+w l) = C ds [ S ds' str \e {1 - s)v We {s - s ' )v We s ' v l\ + ... (34) 
Jo Jo J 

Manipulating this expression further, we will freely use partial integration and shifts in 
/ keeping in mind that (0) appears in the one-loop functional (|30|). Moreover, to avoid 
keeping track of irrelevant terms, we use power counting to argue [cf. Eq. Q2$P ] that the 
divergences bilinear in the fermion fields cannot contain any derivatives of M 2 or fa so 
that we can treat M 2 , fh as constants. Of course, M 2 will not appear at all in the final 
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result. With these qualifications, we obtain for the terms of 0(ipip) 



1 rs 



str (e v+w l) = 2fir 2 g 2 \ ds \ ds' \ e-d^+^X^TM^^-f.-^Mi^t^m)]^ 
- tr 



o Jo 

(l- s+ s'Mi^/+m)^ cie -(s-s')l(l-i^d) 2 +rM'^^ 



2 f ir 2 g 2 J^dsf ds'^ | e -(— ')M(^/+™) e -(i-,+/)[(/+^a) 2 +rM 2 ] 

+e -(l-s+s')^i y ^/+Tm) e -(s-s')[(l+i^d) 2 +rM 2 ]'\^ 



= 2fMT 2 g 2 dz^e-wW^e-V-'W^V + ™ Vd • (35) 
J o 

After integration over z, the one-loop functional bilinear in the external fermion fields is 
found to be 

r r°° dr a f d d l — e -[(i+iV^d) 2 +tM 2 } _ e -^(i^f+Tm) 

Zl=i\w = tm 2 / d x / — r 2 -2 / ^ — 7~=~r~ _> ci • 

^ J Jo t J {2n) a (I 4- lyfrd) + rM 2 — \l\isJt / + rm) 

_ i ( 36 ) 

Except for derivatives of M 2 and m, this is still the complete one-loop functional bilinear 

in the fermion fields. 

To extract the divergent parts, we change integration variables and decompose the 
functional fl36|) in the following way: 

Idl ( ™ J2-i„-TM 2 „ -I 2 



dt t 3-d e -tm e -i/ ) 

ipd • (37) 



o t (i + itd) 2 + t 2 M 2 -/j(it/ + t 2 m) 

The first (/i-dependent) term in this expression is cancelled by the /i-dependent parts 
of the second one. The divergent parts are now easily isolated by expanding the second 
denominator (setting /x = 0) in t. With I := VP, the divergent parts are given by 

~Hiv i 9 f ,h f°° dt n j f d d l — ( 2t sin I _ cos / \ , 



g 



2 



o t i (2tt)^ c1 V 

d 4 x^ cl (- $ + 2m)ip cl . (38) 



(47T) 2 (d-4) 

For the final coefficients we have used the integrals 

^/WsinA = 1 



7 (2vr) d \ I 2 ' I 3 ) " (4 



7T 



Comparing with the divergence Lagrangian (|29|) , we find complete agreement with (|38|) 
in the fermionic part. 
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6. In summary, we have shown that the one-loop functional of a general quantum 
field theory, not necessarily of the supersymmetric type, can be written in terms of the 
Berezinian of a supermatrix operator. We have illustrated - for the case of a Yukawa 
theory - that the divergent parts of the one-loop diagrams may be worked out quite 
easily. The method simplifies in a significant manner calculations in more realistic non- 
renormalizable theories 0. 
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